FERMILAB-PUB-XX-YYY-E, BNL-XXXXX-YYYY-EE, arXiv:hep-ex/XXXX.XXXX 



An improved measurement of muon antineutrino disappearance in MINOS 
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We report an improved measurement of disappearance over a distance of 735 km using the 
MINOS detectors and the Fermilab Main Injector neutrino beam in a X7,j-enhanced configuration. 
From a total exposure of 2.95 x 10^" protons on target, of which 42% have not been previously 
analyzed, we make the most precise measurement of Am^ = [2.621°;^^ (stat.)±0.09(syst.)] xlO"^eV^ 
and constrain the mixing angle sin^(2S) > 0.75 (90%CL). These values are in agreement with 
Am^ and sin^(2^) measured for v^, removing the tension reported in 

PACS numbers: 14.60.Pq, 14.60.Lm, 29.27.-a 



Observations of neutrinos and antincutrinos created in the Sun, the Earth's atmosphere, reactors and acceler- 
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ators provide strong evidence that neutrinos un- 

dergo transitions between their flavor eigenstates {ve, 
Vr) as they propagate. These transitions can occur due 
to quantum mechanical mixing between the neutrino fla- 
vor and mass (yi, 1^2, vz) eigenstates. The mixing may 
be parametrized with a unitary matrix J7pmns [HI which 
is typically expressed in terms of three mixing angles ^12, 
^'23, di3 and a charge-parity (CP) violating phase 5. This 
interpretation, referred to as "neutrino oscillations," re- 
quires that neutrinos have mass and motivates extensions 
to the Standard Model (SM) of particle physics. Exten- 
sions which explain the origin of neutrino masses, for 
example the addition of right handed sterile neutrinos 



12l |. may also explain the baryon asymmetry [13[ of 



the universe. 

The CPT symmetry of the SM requires that 



and 



Vfj^ have the same masses and mixing parameters. In 
vacuum, the probability P{v^ — > v^) that a is detected 
after a distance L as a (rather than a i/g or Vt) must be 
equal to the corresponding probability P{p^ v^j) for 
antineutrinos. For a with energy E the probability 
may be written as 



l-sin2(26')sin2 



1.267ATO2[cV2]L[km] 

E\GcSr\ 



(1) 

where Am^ and s\v?{20) are effective parameters that 
are functions of the angles parameterizing C/pMNS and 
the differences in the squared masses Am?- = — rn] 
of the 1^1, V2 and 1/3 states 
stratcd |A?77|j| 3> 
6*13 « we have 9 
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Experiments have demon- 
14| . In the limiting case that 
and IAto^I ~ sin^(26'i2)| Am|J 4- 
I . Muon antineutrino oscillations 
are described by an equation which has the same form as 
Eq. [T]with parameters Am^ and sin^(26'). The extended 



cos2(26'i2)|Ami2l [15 



SM predicts Am^ = 
vacuum oscillations 



Am-^ and sin^(26i) = sin^(26') for 



mm 



Observation of Pivf^ 



Vfj.) would therefore be evidence for 



physics beyond the SM, such as neutrino interactions in 
the earth's crust that do not conserve Icpton flavor. 

In this Letter we describe a measurement of PiV^ — 5- 
V^) conducted over a baseline L = 735 km using a V^- 
enhanced beam with a peak energy of 3 GeV. The beam 
was produced by directing 120GeV/c protons from the 
Fermilab Main Injector onto a graphite target to produce 
TT / K mesons that decay to produce neutrinos. Two mag- 
netic horns focus the mesons, allowing us to control the 
energy spectrum and I'/D content of the beam. 

The neutrino beam is pointed towards two detectors, 
referred to as Near and Far. The 980 ton Near Detec- 
tor (ND) measures the t'p and content of the beam 
as a function of energy at a distance of 1.04 km from 
the tt/K production target. The 5.4kton Far Detector 
(FD) is located in the Soudan Underground Laboratory, 
734 km from the ND, and remeasures the beam com- 
position. The neutrino detectors are steel-scintillator. 



tracking-sampling calorimeters optimized to identify and 
measure the energy of muon neutrinos and antineutrinos 
and reject backgrounds from neutral current and z/g in- 
teractions [1^. The detectors are magnetized with an 
average field of 1.3 T to distinguish i/^ from 77^ based on 
the charge of the /z produced in weak interactions. 

We previously reported oscillations with an en- 
ergy dependence consistent with Eq. [T] and Am^ ~ 
2.32l2^J^(stat.-^-syst.) xlO-3cV2, sm'^{29) > 0.90 (90% 
CL) 01 . The measurements utilized 7.25 x 10^° protons 
on target (POT) of data collected between 2005-2009 
with a t'^-enhanced _beam (20| . Measurements made by 



Super-Kamiokande Q and T2K 21 1 are in good agree- 
ment with our values. 

In 2009-2010 we collected 1.71 x 10^° POT in a 77^- 
enhanced beam [20j created by reversing the polarity 
of the horns. The magnetic fields in the FD and ND 
were also reversed to focus the created in 77^ interac- 
tions. These antineutrino data also exhibited oscillations 
in agreement with Eq. [U but with parameters Am^ = 
[3.36lJ5;4^(stat.) ± 0.06(syst.)]xl0-3eV2 and sm'^{29) = 
0.86±0.11(stat.)±0.01(syst.) [H. We found no systematic 
effects which could explain the difference between the z/^ 
and parameters. Assuming identical true values for 
(Am2,sin^(20)) and (Am^, sin2(26')) , we calculated that 
such a difference would occur by random chance about 
2% of the time. To clarify the situation we collected an 
additional 1.24 x 10^" POT with the F^-enhanced beam 
during 2010-2011. We have also updated the analysis to 
improve the sensitivity to (Ato^, sin (29)) , reduce uncer- 
tainties due to Monte Carlo (MC) modeling in the ND, 
and increase the similarity to the oscillation analysis. 

We isolate a sample of Vf^^ and charged-current 
(CC) Vfj.N fjiX events by searching for interaction ver- 
tices inside our detectors with a muon track and possi- 
ble hadronic activity from the recoil system X. We re- 
ject hadron tracks reconstructed in neutral-current (NC) 
events by combining four topological variables describ- 
ing track properties into a single discriminant variable, 
/ilD, using a fc-nearest-neighbor (kNN) technique [2^ . 
The kNN algorithm calculates the distance in the four- 
variable space between each measured event and an 
ensemble of simulated events; the output is the frac- 
tion of signal in the fc = 80 closest MC events. This 
discriminant was used in our previous analyses [l|, |^ 
and, as shown in Fig. [Ija), is well modeled by our 
MC. We maximize the statistical sensitivity to Ato^ 
by requiring ^ID > 0.3, which results in a MC esti- 
mated efficiency /purity of 90.7%/99.0% at the ND and 
91.6%/99.0% at the FD. We then discriminate v^, from 
on an event-by-event basis by analyzing the track cur- 
vature in the detector's magnetic field. Figure[T{b) shows 
the track charge/momentum {q/p) divided by its uncer- 
tainty {'7{q/p)), as determined by our track reconstruc- 
tion algorithm. We select j7„ and reject Un by rcquir- 
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FIG. 1: Event selection variables from the ND: (a) The vari- 
able, /ilD, used to select u^/Vfj^-CC events and reject NC 
events, (b) The reconstructed track charge/momentum q/p 
divided by the uncertainty cr{q/p) reported by the track re- 
construction algorithm. In both figures, we have applied 
all selection criteria except the one on the quantity being 
shown. Shaded bands show systematic uncertainties. The 
background histograms are stacked on top of each other. 
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FIG. 2; The V^-CC energy spectrum measured by the ND. 
The selection procedure described in the text has been ap- 
plied. Each point represents one month of data. The shaded 
Monte Carlo band shows the combined effect of systematic un- 
certainties due to cross-sections, flux, energy scale and other 
sources. 



ing ^^gfp^ > with a MC estimated efficiency /purity 
of 98.4%/94.7% in the ND and 98.8%/95.1% in the FD. 
The Vfj^ background accepted by the selection is predom- 
inantly due to high energy muons with small curvature. 

We reconstruct the neutrino energy by summing muon 
and hadronic shower energies. The muon energy is mea- 
sured using track range and curvature. We reconstruct 
the hadronic shower energy using three variables: the 
sum of the reconstructed energy deposited by showers 
that start within 1 m of the track vertex; the sum of 
the energy in the two largest showers reconstructed in 
the event; and the length of the longest shower. We use 
these three variables in a second kNN algorithm and es- 
timate the shower energy as the mean true hadronic en- 
ergy of the k = 400 closest MC events. This technique 
improves the hadronic energy resolution when compared 
to a method which uses only the energy deposited by 
the largest shower, increases the statistical sensitivity to 
by 10%, and was previously used to analyze dis- 
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appearance 

Data from the ND are used to predict the neutrino en- 
ergy distribution at the FD. Though both detectors have 
the same segmentation and very similar average magnetic 
fields, for economic reasons the ND is smaller and asym- 
metric about the magnetic field coil and is more coarsely 
instrumented with scintillator in the downstream "muon 
spectrometer" region [l^ . In addition, the ND coil occu- 
pies a larger fractional area than the FD coil and more 
muons enter it. In the ND data, we observe a reconstruc- 
tion failure rate of 6.1%, mostly associated with tracks 
entering the coil region, but the MC predicts 4.2%. Pre- 
viously, we dealt with this issue by assigning a systematic 
error. Now, we remove ND events with a track that ends 
less than 60 cm from the coil. We also remove events 
with a track that ends on the side of the coil opposite 
the beam centroid. The new event selection decreases the 
efficiency to 53% in the ND, but reduces the data/MC 
failure rates to 1.4%/0.9%. The selected sample contains 
the same classes of neutrino scattering processes as are 
present at the FD, and our results are not significantly 
more vulnerable to cross-section uncertainties. We ap- 
plied the new selection and shower energy reconstruction 
to the 2009-2010 data and found that the best fit pa- 
rameters shifted by only (5(Am^) = +1.0 x 10"'* eV^ and 



(5(sin^(26l)) 



-3.6 X IQ- 



A total of 2.98 x 10^° protons were delivered to the 
graphite target during the data-taking periods. We 
impose a series of data and beam quality 0, [ij, [IH 
requirements which reduce the analyzable exposure to 
2.95 X 10^0 POT (99.0% livetime) at the FD, and 2.73 x 
10^0 POT (91.8% livetime) at the ND. The uncertainty 
in the livetime is negligible, and the high and largely 
overlapping livetime in both detectors assures that our 
results are not sensitive to beam effects that would cause 
the number of events per POT to vary. 

Under normal conditions, the ND measures about 2400 
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FIG. 3: The reconstructed energy of FD 77^-CC events. These 
data are fit to produce the 2009-20if contour shown in Fig. [l] 
The band displays the effect of systematic uncertainties. 



Vfj^-CC events per day in the oscillation energy region 
{El, < 6 GeV). These data are essential for monitoring 
the neutrino beam and the quality of the experiment. 
Figure [2] shows the reconstructed energy distribution 
measured in each month during the data-taking periods. 
With the exception of Feb. 2011 (the last point in each 
bin in Fig. [2]), all months are in statistical agreement, 
and we expect a constant counting rate per POT at the 
FD. Part of the February dataset was taken after the 
neutrino target's cooling system had failed, leaking water 
into the target canister. This resulted in a decrease in the 
neutrino flux of 4% from 0-6 GeV when integrated over 
the entire month. The decrease is adequately modeled 
by our beam simulation, and we account for it at the FD 
using the ND data. 

Measurements of the beam position and width at the 
target, the position of the remnant proton beam at the 
end of the decay pipe, and the integrated muon flux from 
meson decays in the pipe, all further indicate that the 
expected number of F^-CC per POT is constant at the 
FD. Using the proton beam and monitors we measured 
the target and horn misalignments as well as the residual 
magnetic field in the neck of the focusing horns [2^ . We 
conclude that these effects introduce <C 1% uncertainty 
in the rate at the FD. 

Events at the FD are read out in a 100 /zs window 
surrounding the 10 ^s long beam spill 27|. We select 
events coincident with the beam spill that have a ver- 
tex in the fiducial volume and a track identified as a 
muon by ^ID. We remove cosmic-ray background by re- 
quiring that the cosine of the angle between the track 
and beam direction be > 0.6. A total of 521 events sat- 
isfy our criteria, 273 in the 2009-2010 dataset and 248 in 
2010-2011. Using the muon charge, we identify 328 as 



readout windows taken in anticoincidence with the beam 
during 2005-2011 and estimate a cosmic-ray background 
of 0.8 events. We monitor detector and reconstruction 
performance by measuring the rate of cosmic-ray muons 
traversing the detector before and after the beam spill. 
The rates in each data-taking period, 0.38 ± 0.03 Hz in 
2009-2010 and 0.41±0.04 Hz in 2010-2011, are consistent. 

We predict the energy spectrum at the FD by first 
correcting the ND spectrum for inefficiency and back- 
grounds. We then transfer that spectrum to the FD us- 
ing a two dimensional "beam matrix" Q. We predict 273 
events if F^ do not oscillate, including 3.5 NC and 10.8 



CC [31|. The energy spectra are shown in Fig. [31 Os- 



cillations are incorporated into the prediction according 
to Eq. [TJ Maximizing the binned log-likelihood yields 



[2.62i°:^^(stat.) ±0.09(syst.) 
sin^(25) = 0.95ig;i'i'(stat.) ±0.01(syst.) 
sin2(25) > 0.75 (90% C.L.) 



,2^ 



xlO^^cV^ 



j/^ and 193 as i^^. We apply the selection to 37.7 x 10^ 



with p = 35.3% at best fit [32|. When analyzing 
the 2010-11 alone we obtain Am^ = [2.26to:29(stat.) ± 
0.09(syst.)] xl0-3eV2, 3111^(26) > 0.79 (90% C.L.) with 
p = 14.5%. 

The systematic errors stated above were evaluated by 
computing the standard deviation in the best fit param- 
eters as the fits were repeated using MC samples that 
were shifted in accordance with uncertainties in neutrino 
cross-sections; the beam matrix and relative FD normal- 
ization; NC and ly/i-CC backgrounds; the relative FD to 
ND energy calibration; the absolute muon energy scale; 
and the absolute hadronic energy scale, including final 
state hadronic interactions. The input uncertainties are 
as in our previous analyses [l|, |3| and were derived from 
in-situ data, bench tests of detector and beam compo- 
nents, a test beam experiment, and published neutrino 
and hadron cross-sections 0, [l^ . We increased the un- 
certainty in the axial-vector mass from 15% to 30% to 
account for additional uncertainties in F^ quasi-elastic 
scattering pgj . 

w 

The confidence regions shown in Fi g. HI were calculated 
according to the unified procedure of [30| and incorporate 
both statistical and systematic uncertainties. We assess 
the consistency of the and F^j measurements by do- 
ing a joint (Am^ = Am'^ , sin^ {26) = sin^(26')) fit to the 
data to establish the null hypothesis. We then perform 
four-parameter (Ar?!^, Am^, sin^(20), sin^(26')) fits on an 
ensemble of z^^ and F^ MC experiments generated at the 
joint best fit. The joint fit had a larger likelihood than 
p = 42% of the four parameter fits, indicating consistency 
between i^^ and F^ . 

In conclusion, we have used a F^-enhanced Fermilab 
accelerator beam and detectors that discriminate u^, from 
F^ to make the most precise measurement of Arn'^. Our 
results remove the tension reported in [l| and establish 
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FIG. 4: Confidence regions calculated by fitting the FD data 
with an energy spectrum predicted from the ND but modi- 
fied to incorporate oscillations via Eq. [T] The regions include 
statistical and systematic uncertainties. The 2009-2011 
contour is derived from the energy spectrum shown in Fig. [3] 



consistency between and i^fj, oscillations at L/E w 
200km/GeV. 
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